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Abstract — In this research, weighted total acceleration for a function f (x, z, t)was formulated. This total 
acceleration equation was done at the Euler momentum equation. Then, the Euler momentum equation was done 
together with free surface boundary condition equation to formulate water wave constant at the solution of 
Laplace equation. The velocity potential of the solution of Laplace equation actually consists of two components 
that were used in this research. 
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I. INTRODUCTION 

Momentum equation is an important basic equation in 
mathematic modeling of hydrodynamics, including water 
wave modeling. Momentum equation commonly used in 
water wave modeling is Euler momentum equation. There 
is a constraint in this equation, i.e. Euler momentum 
equation has no hydrodynamic force in the horizontal 
direction or convective acceleration has a value of zero 
when velocity potential is substituted to the term. To 
overcome this problem, weighted total acceleration 
equation was formulated where there are two weighted 
coefficients, i.e. at the timet differential term and at the 
differential term of vertical-zdirection. 

Laplace equation solution consists of two velocity 
potential components (Dean (1991)). However, only one 
component that has been used. Equations from water 
wave constant, i.e. wave number fcand wave constant G 
can be formulated using only one velocity potential 
component, but the value is determined by both the two 
velocity components. In this research, the water wave 
surface equation is formulated using the two velocity 
potential components, then the condition of the water 
wave surface that has been produced is studied. 

H. WEIGHTED TOTAL ACCELERATION 

Hutahaean (2019a) formulated weighted total acceleration 
in a function / = f(x,t), x\s horizontal axis and t is 
time, using Taylor series. The formulation of weighted 
total acceleration in a function/ = f{x,z, t),zis vertical 
axis, is done using similar method, therefore the 
formulation of weighting total acceleration in / = 
fix,z, t)will be preceded by reviewing the formulation of 


weighting total acceleration in f — f(x,t) to obtain a 
clearer description. 


2.1. Weighted Total Acceleration for the function off = 
f(x,t) 

The changes in the value of a function in a function / = 
f(x,t )for a very smallbxandbtusing Taylor series only 
until the second derivative is, 

a / a / 

fix + Sx, t + St ) = fix, t) + Sx — —f 3t — 

dx at 

Sx 2 a 2 / a 2 / st 2 a 2 / 

H r—~—x + StSx -1 ————— 

2 3x 2 ata* 2 at 2 

By working on the argument of Courant (1928) that in 
order to obtain a good result on horizontal velocity u = 

dx 

—, then weighting coefficient y, is done which is a 
positive number, in time differential in Taylor series. 

a / a / 

f(x + Sx, t + ySt ) = fix, t) + Sx- —h y3t — 

dx at 

. Sx 2 a 2 f „ „ a 2 / y 2 st 2 a 2 / 

-i-v + yStSx —— +-. . .(l) 

2 ax 2 ' atax 2 at 2 v ' 

At the limit Sx, btclose to zero the following equation is 

obtained, 

.(2) 

This equation is weighted total derivative equation or 
weighted total acceleration for the functionof/ = 
fix, tj where yis weighting coefficient. 


Df a/ a/ Df a/ a/ 

— — u —h y —or— — y —I - u — 
dt cix a t dt at ax 


The method of calculating weighting coeffecient ywill be 
formulated using Taylor series (1). The second derivative 
term can be omitted if, 

Sx 2 a 2 f ... a 2 / y 2 «t 2 a 2 / 

~r^ + y stSx 


Sx ¥c + y at I 


< £ 


•(3) 
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Then it was defined Sx — CSt = - St = — St = - St, 

T kT k 

whereCis wave celerity, k is wave number T is wave 
period, a = ~ is angular frequency. 8x in (3) is 
substituted with ^5t, and the following equation is 


obtained, 


a 2 Sta 2 f , yaSt a 2 / , y 2 Sta 2 f 
2 k 2 ax 2+ k stax + 2 at 2 


era / a/ 


< £ 


■(4) 


The completions of this equation requires a function form 
of f — f(x,t). And the following sinusoidal function 
form will be used. 


f(x, t) = cos kx cos at ....(5) 

This equation is water wave surface equation of the linear 
wave theory. The derivative of the function is as follows 


Table. 1: Derivative Equation of (5) 


a/ 

dx 

= —k sinkx cosat 

a 2 / 

ax 2 

= — k 2 cos kx cos at 

a / 
at 

= —a coskx sin at 

a 2 / 

atax 

= ka sin kx sin at 


a 2 / 

at 2 

= —a 2 cos kx cos at 


Using the condition of coskx = sinkx — cosat = sinat, 
the elements of sinusoidal function will cancel out each 
other as a result of a division. Substitute the derivative 
equations to (4), the following equation is obtained 

— Y + ~Y 

2 2 


1+7 


a St 


The numerator (1 + y)is a positive number, then the 
equation can be written as, 


2~ r + 2 r ‘ 




If equals (=) relation is used, then 

;-)' + ;)' 2 = io + rt .( 5 ) 


Considering that yis a positive number, the right side of 
the equation is a positive number. Therefore, the left side 
of the equation is also a positive number. The calculation 
of the value yean be done by releasing the sign | | in the 

left side of the equation, i.e. using equation (5). 


The calculation of the valueywith (5) requires an input St. 

The value of St. is obtained from the function / = f(t). 

The approximation of Taylor series for the function is, 

a/ 5t 2 3 2 / 

/(t + St)-/(t) + St—+ —— 


In order to be able to be used only until the first 
derivative, then 


st 2 a 2 f 


sta 2 f 

2 at 2 


2 at 2 

St if 


a / 

at at 


at 


< £ 


For the 


function, fit) = cos at: — = —a sinat: —7 = 

’ 1 v J ’ at ’at 2 

—a 2 cos otand it is done in a cosat — sinat condition, 

and 

Sts 2A 

--< e or ,ot = —.( 6 ) 


is obtained . Substitution of ( 6 ) to (5) obtains 
7 = 3 ...(1) 

It is obtained that yhas a constant value, i.e. independent 
of wave period or the level of accuracy e. 


2.2. Weighted Total Acceleration for the function / = 
fix, z, t) 

To obtain weighted total acceleration equation in a 
function / = fix, z, t), the similar method will be done as 
in the function / = f(x,t ), where, 
fix + Sx, z + y z Sz, t + y5t) = 

a / a / a / 

fix, t) + Sx — + y z Sz — + y3t — 

' ax az at 

J x2 z 2 f^ xx a 2 / , (7 z 5z) 2 a 2 / 

H —~~—— + y z SzSx —— + 


3z 2 


+Y S tS x^+YY z St8z^ + ^ a 4 .. 

' ata* ,z ataz 2 at 2 

In ( 8 ), foiAz = Sx, it is meant that y z Sz^- — Sz (y z 

therefore in a change of z for Sz = Sx, the value of the 

first derivative function against zis (y z Jy), and so also 

y z 8z8x —^-is meant5z5x fy z —^-)and (yzl?z) 3 { which 

,z azax V' z azax/ 2 az 2 

Sz 2 ( 7 a 2 /\ * • 1 • -i 

means as 5^2 J- m th e previous section, the 

value of5xand5zis5x — Sz — CSt — - St — = —. 

T kT k 

Then, a function f(x,z,t) is reviewed with the following 
form. 

f(x,z,t ) = coskxcoshkifi + z) cosat....(9) 

Atz = 0, c 1 = cosh(kh)andc 2 = sinh(kh) are defined 
and done in the deep water where tankh — lwith the 
value of kh — 2.0n. Thenc x = cosh(2.0n ) — c 2 = 
sinh(27r),and ( 8 ) is done in a condition of coskx = 
sinkx — cosat — sinat, then the sinusoidal function 
cancelled out each other. The derivative equations (9) can 
be written in the forms shown in Table (2). 

Table. 2: Differential of (9). 



= — k 2 Ci 
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az 2 

a 2 / 

3x3 z 
= -k 2 c 2 

a 2 / 

3t3z 
= — okc 2 

a / 
at 

= —oc 1 

32/ - k 2 c 
az 2 “ k Cl 

a 2 / 

at 2 = _fJ Cl 


To simplify the writing, the followings are defined 

Sx 2 a 2 / a 2 / ( y z Sz ) 2 a 2 / 

A = ^d + ^ 5 zSx zd- x + ^d 

a 2 / a 2 / oat) 2 a 2 / 

+yStSx ———h yy z StSz „ - + ■ 


3t3x 


3t3z 


at 2 


a/ a/ a/ 

B - Sx — + y z Sz — + y3t — 

ax /z az r at 

In order for ( 8 ) to be able to be used with only the first 
derivate, then 

1-1 < £ .(9) 

IbI 

The substitution of differential equations in Table (2) to 
(9) will obtain, 


y , 

r* +r 


—-(- 

2 aSt K 


■yc i - cO 


-( r+ 1 + ^) c ^ + i r ’ =0 

2e 

Substituted'tfrom ( 6 ), St — — 


y , 
'Y Cl + r 


2 2 


-yci - c x ) 

- (y + 1 + ^) c 2 y z +jYz = 0.(10) 

With (10), y z can be calculated whereyis a known 
from(7). With an input y = 3, y z — 1,630 is obtained for 
Ci = cosh(2.07r)andc 2 = sinh(2.0n ) where c 1 = c 2 . 

As a result the second derivative in ( 8 ) can be omitted and 
the total derivative equation for function / = /(x,z, t)is 

— -y — +u— + y 7 w— .(11) 

dt 'at a* ,z az v ' 


III. A Complete Velocity Potential Equation 
By completing the Laplace equation with separation 
variable method and after doing the time periodic 
boundary condition and lateral periodic boundary 
condition. Dean (1991) obtained velocity potential 
equation that consisted of two potential velocities, i.e. 

<p(x, z,t) — A cos k x(Ce kz + De ~ kz ) sin a t 
+B sink x(Ce kz + De~ kz ) sin a t.. .(12) 
cp{x,z,t ) = <p A (x,z,t ) + (p B (x,z,t)...( 13) 

cp A {x,z, t) = A cos kx(Ce kz + De~ kz ) sin a t 
...(14) 

cp B {x,z, t) = B sin k x(C.e kz + De ~ kz ) sin a t 
...(15) 

There are four constants that should be determined, i.e. 
A, B, CandD. Hutahaean (2019b) has shown that the two 
equations have similar constant value, or in other words 


there is only one constant value in velocity potential total 
(12). However, in the next section it will be proven again 
with another method that ( 12 ) has one constant value. 

Equation (12) can be written as, 

cp{x,z, t ) = (A coskx + B sinkx ) 

( Ce kz + De~ kz ) sin a t .(13) 

At a condition of coskx — sinkx, (13) can be written as 
cp(x,z, t) — (A + B ) cos k x 

( Ce kz + De~ kz ) sin a t .(14) 

or 

cp(x,z,t ) = (A + B) sinkx 
(Ce kz + De~ kz ) sin a t .(15) 

The constants of A, B, CandDwill be formulated using 
(14) and (15), where it will be proven that either using 
(14) or (15) similar constant will be obtained. The 
formulation is done by doing kinematic bottom boundary 
condition on flat bottom, as was done by Dean (1991). 

a. Alternative I 

The constants^, B, CandD will be determined using (14) 
where water particle velocity at the vertical-zdirection is 
dtp 

w = — — = —(A + B)k coskx 
oz 

( Ce kz — De~ kz ) sin a t 

Substitute equation for w to the kinematic bottom 
boundary condition equationw_ h = —u_ h where at flat 

bottom — = 0 , 

ax 

— (A + B)k cos kx(Ce~ kh — De kh ) sin at = 0 

The equation is divided by — (A + B)kcoskxsinatfor 

coskx =4 Oand sinat =4 0 

Ce~ kh — De kh — OorC = De 2kh . Substitute C (14) 

<P(x, z, t) = (A + B ) cos k x 

(De 2kh e kz + De~ kz ) sin a t 
or 

0(x, z, t) = (A + B)De kh cos k x 

( e fc(h+z) + e - m+z) )sina t 
A new constant is defined 

G a = (A + B)De kh .(16) 

0(x,z, t) = G a cos k x{e k(jl+z ^ + e -k(h+z) ) sin a t 

.(17) 

a. Alternative II 

The constants A, B, CandDwill be determined using (15), 
dtp 

w — — — = — (A + B)k sin k x 
dz 

( Ce kz — De~ kz ) sin a t 

Substitute equations for u and w to the kinematic bottom 

boundary condition equation 

dh , ah - 

w_ h — —u_ h —, where — =0 

n n dx ax 
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— (A + B)k sin k x(Ce~ kh — De kh ) sin at — 0 

The equation is divided by — (A + 

B)ksinkxsinatfoxsinkx & Oand sinat & 0 

Ce~ kh — De kh — OorC = De 2kh . Substitute Cto (15) 

0(x, z, t) = (A + B) sin k x 

(. De 2kh e kz + De ~ kz ) sin a t 
or 

0(x, z, t) = (A + B)De kh sin k x(^e k ^ h+z ^ 

+ e~ k ( h+z ^ sin a t 
A new constant is defined 

G b = (A + B)De kh .(18) 

<P(x, z, t) = G b cos k x 
(e fc(ft+z) + e -W+z^ sin a t 19) 

From(16) and (18) obtained that G A — G B — G , so it is 
proven that in (1) there is only one wave constant 
valueG.then (7) becomes 

0 ( x, z,t) — G (cos k x + sinkx ) 

^gk(ft+z) _|_ g-fc(h+z)^ s ^ n a j ( 20 ) 

The hyperbolic function equation is, e klh+z> + 
g-fe(/i+z) _ 2 coshkfi + z), (13) becomes 
0(x,z, t) = 2G(cos kx + sinkx)coshk(h + z) sin a t 
Defined G — 2G 

0(x,z,t) = G ( cos kx + sinkx) coshk (h + z) sin a t 
....(21) 

A complete velocity potential equation is obtained with 
the form as in (21). In that equation, there are still two 
wave constants where the form should be known, i.e. 
wave number /rand wave constant G. Considering that 
the values of wave number /rand wave constant Gis 
similar along the wave curve, then the calculation of the 
two parameters will be done at the point of characteristic 
where coskx — sinkx, at this conditional) becomes, 
0(x, z, t) = 2Gcoskxcoshk(h + z) sin a t....(26) 


The particle velocity in horizontal-xdirection is, 
u — — ^ = 2Gksinkxcoshk(h + z)sinat ....(27) 
The particle velocity in vertical-zdirection is, 
w = — — = —2 Gkcos kx sinh k(h + z) sin at .(28) 


IV. Application of Weighted Total Acceleration on 
Euler Momentum Equation 

From (28), the total derivative for horizontal ^direction 
velocity is, 

Du 3 u 3 u 3u 


-T- = y — + U -—b v 7 w — 
dt ' 3t 5x 3z 

With this total derivative equation, the Euler momentum 
equation in horizontal-^direction becomes, 

3 u 3 u 3 u 13p 

V~ r+ u — + y z w— = — 

3t dx 3 z pdx 


By doing the characteristic of irrotational flow. 


3u 

dz 


3w , . t 

—obtained, 

ax 

au i a , 2 , 2 \ i ap 

y - 1 -(tr + y z w 2 ) =- 

'at 2 ax K ,z J pax 

Total derivative equation for vertical velocity in axis-z 

direction. 


.(29) 


Dw 


3w 


3w 


3w 


W = r at +u a^ +r - u ’aF 

The Euler momentum equation in vertical-zdirection 
becomes, 

3w 3 w 3w 13p 

p 3 z 


r IF +u ^ +r - w H 


9 


The execution of irrotational flow characteristic, — = — 

a% 3z 


3w , 1 3 r 2 i 2 a 

y ir + + ^ ) = 


1 3p 
p dz 


g .(30) 


(29) and (30) are modified Euler momentum equations, 
where there are time weighting coefficient yaiul 
weighting coefficient vertical zdirection of weighting 
coefficient, i.e.y z . Using (30) pressurep equation will be 
formulated where (30) is written as an equation for 
pressurep. 

13p aw i a 

-W^ y 3t + ~23i iU+V ^ ) + B 

This equation is multiplied by dz and integrated against 
vertical-z axis. 

P = 

p r J z at ~ ' 2 

1 

- r (u 2 + y z w 2 ) + g(r / - z) 


? ' aw !/, « 
dz + - (u 2 + y z w 2 ) 


Differentiated against horizontal-xaxis 
13p 3 f"3w 13 , , 

W* = y ail lt iz + 2a;^- + '-“'O 

1 a 3 tj 

(U + y z w 2 ) + g 


2 ax 

Substituted to (29) 

3u 

y 


ax 


3u 3 flaw 

at + Y a^) z aF dz 


dx 

The completion of — f'J dz\s done using velocity 

potential (21), where the particle velocity in horizontal 
direction is in equation (27), and the particle velocity in 
vertical-zdirection (28). From (28) the following is 
obtained, 

3w 

= —2Gkacoskxsinhk(h + z)cosat 
This equation is integrated against time t. 


I 


11 3 w 

——dz = —2Gacoskx 

at 


(i coshkfi + p) — coshk(h + z))cosat 
Then, it is differentiated against horizontal-* axis 
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a f'aw 

— —dz = Gkasinkx 

3 xj z at 

(i coshkQi + ?;) — coshk(h + z))cosat 

Equation (27) is differentiated against time t, — = 

at 

2 Gkasinkxcoshk(h + z) cosat, and it is seen that this 

form is in — [ v — dz, so the following relation is 
ax J z at b 

obtained 
— [ ? — dz = 

3x J z at V at at/ 

Substitute this equation to (31) 


^17+ “K+ Y z w*) = -g% .(32) 


This equation is a surface momentum equation that will 
be used in the calculation of Gandk. 


V. THE FORMULATION OF AN EQUATION 
FOR THE CALCULATION OF G AND k 

As has been mentioned in the previous section that the 
calculation of Gand kis done in the point of characteristic 
where coskx — sinkx. Therefore, (27) is used as the 
particle velocity in horizontal ^direction and (28) is 
particle velocity equation in vertical z direction. 


5.1 Wave number conservation equation 
In the formulation of an equation for the calculation of 
Gand/cin the following sub-chapter, the wave number 
conservation equation will be done. The equation come 
from the principle of variable separation at the completion 
of Laplace equation, i.e. that velocity potential is 
considered as a multiplication of three functions, i.e. 
<P(x,z,t ) = /C(x)Z(z)T(t)whereA'(x)is just a function-x, 
Z(z)is just a function-zandT(t) is just a function-t. In 
this case Z(z) = coshkQi + z). As just function-zthen, 

3 Z(z) 


3 coshk(h + z) 3 k(h + z) 

-= sinhk(h + z)-= 0 

3x dx 

For sinhk(h + z)is not equal to zero, then 
^1^=0 ....(33) 


This equation (33) is called wave number conservation 
equation. This means that all area of calculation has 
similar values for the function tanhk(h + z), coshk(h + 
z) dan sinhk(h + z). As deep water, it can be defined as 
water depth wher etanhk(h + ?;) = 1 , where?; = ?;(x. t)is 
the water surface elevation against still water level. 
Bearing in mind that the wave number conservation 
equation or law, in the entire domain applies 
tanhk(h + rf) = 1 ....(34) 

In this research, the following is used 
k(h + rj) - 2.0 n ....(35) 

Wheretan/i(2.07r) = 0.999993. 


sinhk(h + 77 ) = coshk(h + 77 ) = sinh(2.0n ) = 
cosh(2.0n ).(36) 


Bearing in mind this wave number conservation law, then 
even though the weighting coefficient is formulated in 
deep water condition, it will also apply in other depths. 

5.2. Substitute Velocity Potential to Momentum Equation 
From (27), 

u = 2Gksinkxcoshk(h + z)sinat 
3u 

— = 2 Gk coskxcoshkih + z)sinat 
dx 

3 u 

u — = 4 G k sinkxcoskxcosh k(h + z)sin at 
dx 

At z = ?; 

u — — 4G 2 k 3 sinkxcoskxcosh 2 k(h + ri)sin 2 at (37) 


From (28) 

w = —2Gkcos kx sinh k(h + z) sin at 
3 w 

—— = 2 Gk 2 sin kx sinh k(h + z) sinot 
3x 

a w 

w-^— = —4 G 2 k 3 sinkxcoskxsinh 2 k(h + z)sin 2 at 
At z = ?; 

w— = —4 G 2 k 3 sinkxcoskxsinh 2 k(h + rj)sin 2 at) 


..(38) 

In deep water wher etanhk(h + 77) = 1, then sinhk(h + 
rj) — coshk(h + r]). Substitute (37) and (38) to the 
convective velocity, 

/ 3 u 3w\ 

("S + 7 '“'e7„„ = 

(1 — y z )G 2 k 3 sin kxcoskx cosh 2 k(h + 77) sin 2 at 
At the characteristic point, i.e. a point where coskx = 

A 

sinkx — cosat — sinat, where?; = -, 

/ 3 u 3w\ 

( u ^ +wr ^)„„ = 

^ (1 — Yz)G 2 k 3 cosh 2 k (ji + 0 .(39) 

If in (39) y z = 1 is used, then ~ + w Yz = 0 

z —?7 

will be obtained. So, it is found that if at the term 
w —andy z = lis used or without weighting coefficient y z 


, hydroynamic force of the surface in horizontal 
direction has a value of zero or there is no 
hydrodynamicforce. From (27), at the characteristic point, 

^ = Gka cosh k(h+^j .(40) 

Substitute (39) and (40) to (32), 

yaGk cosh k (h + —^ 

+ \ (1 - Yz)G 2 k 3 cosh 2 k (h + 0 = -g ^ ....(41) 

Where g ^-is worked at the characteristic point. This 
equation is a relation between Gand/cwhere wave 
amplitude Ais the input. 
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5.2 The formulation of wave amplitude function 

The weighted total acceleration equation (2), done at the 

water wave surface equation p — rj(x,t), obtained^-= 

y — + u 1} The original kinematic free surface 

' at ^ ax & 


a?; 


boundary condition (KFSBC) equation is, w r] = + 


ar; 


u v —. By comparing the two equations, then the KFSBC 


1 ax 

equation should be in the form of = y ^ + u v ^ , or 

y — = w„ - u„ — ...(42) 

'at v ’'ax v ’ 

Substituteufrom (27) andwfrom (28) and done at z = rj. 

dp 

y — = —2Gksinhk(h + r/)coskxsinat 
3 1 

—2Gkcoshk(h + r])sinkxsinat ^ .(43) 

Water wave surface equation was obtained by integrating 
(43) against time t. The right side of the equation is a 
non-linear function against time tof which the integration 
completion is difficult. However, there is an argument 
that can simplify the integration (43) completion. First 
bearing in mind (36), i.e. coskh(h + rf) = cosh(2.07r) = 
constant. Then, (43) is written as, 
dp 

Y —— -2 Gk 
r dt 


3m . 

) — sinat 


/ dm 

I coskxsinhkQi + rf) + sinkxcoshkQi + rf) — I : 

...(44) 

In (44) the one that is the function of time tis only the 
element of sinot. In addition, as a periodical function 
against time t, the element —2Gk [coskxsinhkQi+ 

p) + 2sinkxcoshk(h + 77 )^ should be a constant 

number against time t. Thus, the integration (44) against 
time t, is sufficient by integrating only the sinertelement, 
obtained 


2 Gk 


( dm 

[coskxsinhkQi + p) + sinkxcoshkQi + p) — I cosat 


.(45) 


At the characteristics point, (45) can be written as 

2 Gk 

( dm 

\sinhk(h + p) + coshkQi + p) — I coskxcosat 

.(46) 

The form coskx was selected becauseit has been 

determined that the velocity potential component that was 

used is coskx componentlt is defined 

2 Gk ( dm 

A — - sinhkih + n) + coshk(h + rj ) — 

ya V dx/ 

Then (46) becomes 


r/(x,t) = Acoskxcosat .(47) 

At the characteristic point, then)/ = -, wave amplitude 

function equation, 

2 Gk f 


2Gk ( ( A\ I A\kA\ 

A =- sinhk [h -1 — — coshk [h -1 — — 

ya V \ 2 / \ 2)2) 

From (36) wherestnh/c + ~) = coshk (h + 0 the 

wave amplitude function equation becomes, 

A — — coshk 


ya 


(* + 0(l-“).W 


5.3 Equation for the calculation of /cand G 
Substitute (47) to (41) at the characteristic point 

yaGk cosh k (h + —^ 

+ ^(1 - y z )G 2 k 3 cosh 2 k (h. + = g 

Substitute wave amplitude function, 

/ A\ 1 

Gya cosh k I h + — J + — (1 — y z )G 2 k 


kA 

~2 


cosh 2 k 



Gk 

( A\f 

^ + 2 ) : 

= a —cosh k 
ya 



kA\ 

YJ 


The equation is divided by — cosh k (ji + 

(1 — y z )Gk 2 cosh k (h + — ^ = gk ^ 


2 2 Y a 

ya +y 


kA\ 


Wave amplitude equation is written as an equation for G , 
i.e. 

Ayer 


G = 


-(49) 


2kcosh k ( ,l +f)( 1 -^) 

and substitute it to the last equation. 


Y 


2 a 2 ( 
.(50) 


kA\ y 2 a 2 

1 ~y) + ~y 


(1 - y z )kA = gk(^ 


kA\ 


The calculation of the value A:with this equation using 
Newton-Rhapson method requires initial estimation 
of/efor the initial value of the iteration. The initial value of 
/ccan be obtainedby ignoring convective acceleration, then 
(50) becomes 

y 2 a 2 =gk(l- y) ....(51) 

This equation is the quadratic equation of wave number 
/cthat can be easily completed. The use of (51) maximum 
value of wave amplitude A in a wave period in deep water 
is obtained, i.e. if the determinant D from (51) has a value 
of zero. 


A — 


2y z (7 z 


....(52) 


The value ofGcan be calculated using (49). 

VI. THE FORMULATION OF WATER WAVE 
SURFACE EQUATION. 

Water wave surface equation is formulated using a 
complete velocity potential equation, i.e. equation (21). 
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By using (21), particle velocity in horizontal-indirection 
and particle velocity in vertical-zdirection are 
consecutively, 

u = Gk(sin kx — coskx)coshk(h + z) sin a t 
w = —Gk(cos kx + sinkx)sinhk(h + z) sin a t 
The two particle velocity equations are done at z = pand 
substituted to equation KFSBC (42), 
dp 

Y — = —Gk(cos kx + sinkx)sinhk(h + p) sin a t 
at 

ap 

—Gk(sin kx — coskx)coshk(h + p) sin at — 

dm 

....(53) 

As in the previous section, the water wave surface 
equation is obtained by integrating (53) against timet, 
where the integration is sufficient to be done only at the 
sin a telement. 


j?0,t) = 

Gk 

— (cos kx + sinkx) sinhklh + p) cos a t 
Y a 

Gk dp 

H-(sin kx — coskx)coshk(h + p) cos at — 

per 3x 

In the deep water the equation can be written as, 

V(x, t) = c 0 ((c 2 + Cj) + (cj - c 2 ) ^ c 3 .(54) 

where.to simplify the writing c 0 = — coshk(h + p), c 1 = 
sinkx, c 2 — coskx dan c 3 = cosat are defined. Equation 


(54) is differentiated against horizontal-* axis 
^ = c 0 k ((-Ci + Ci) + (c 2 + cO ^ c 3 .(55) 


Equation (54) is water wave surface equation that is used 

to calculate water surface elevation where —in (54) is 

ax ' 

calculated using (55). pin c 0 ——coshk(h + rj) is 

calculated using the equation, 

p(*, t) = A(coskx + sinkx)cosat ....(56) 

Whereas —in (55) it is calculated with, 

— — Ak(—sinkx + coskx)cosat ....(57) 


Table. 3: The result of calculation of wave parameter and 
other characteristic 


T 

H 

L 

H 

H 

V max 

(sec) 

(m) 

(m) 

1 

A 

H 

6 

1,409 

5,026 

0,28 

2,865 

0,851 

7 

1,918 

6,842 

0,28 

2,865 

0,851 

8 

2,506 

8,936 

0,28 

2,865 

0,851 

9 

3,171 

11,309 

0,28 

2,865 

0,851 

10 

3,915 

13,962 

0,28 

2,865 

0,851 

11 

4,737 

16,894 

0,28 

2,865 

0,851 

12 

5,638 

20,105 

0,28 

2,865 

0,851 

13 

6,617 

23,595 

0,28 

2,865 

0,851 

14 

7,674 

27,365 

0,28 

2,865 

0,851 

15 

8,81 

31,413 

0,28 

2,865 

0,851 


Using water wave surface equation, the elevation of wave 
crest p ma xand the elevation of wave trough p m i n are 
calculated. The wave height is H — ri max — ri min , 
whereas Wilson (1963) criteria i S ?/ '" ax .'fable (3) presented 
the result of the calculations of wave height, wavelength, 
wave steepness, and the comparison of wave height H and 
wave amplitude A. 


Wave-steepness ^ = 0.280, where considering the 
calculation used maximum wave Amplitude A that was 
calculated using (52), then wave steepness is critical wave 
steepness. 

Table.4: Types of wave, according to Wilson criteria 


(1963) 


Wave Type 

V max 


H 

Airy waves 

< 0.505 

Stoke’s waves 

< 635 

Cnoidal waves 

0.635 < T]max < 1 

H 

Solitary waves 

= i 


VII. THE RESULTS OF THE EQUATION. 

7.1 The characteristic of water wave surface. 

In the calculations that will be done in this section, the 
value of y — 3.0andp z = 1.630are used and the 
calculation is done in the deep water. Deep water depth 
h 0 is obtained with the following equation 

h 0 (2.0tt-“) ....(58) 

Where/1 is calculated using (52). 


The critical wave steepness is bigger than the criteria 
ofMichell (1893) i.e.j = 0.142. The comparison 

between wave height and wave amplitude isj = 
2.865which is bigger than 2. Therefore, therelation 
between wave height and wave amplitude is H — 2A 
cannot be used. The obtained Wilson parameter is Vm “ x — 

0.851. Based on Wilson criteria (1963), Table (4), the 
value of the parameter shows that the wave profile has a 
cnoidal wave type, with wave profile presented 
inFig.l.and Fig.2.for wave period T = 8 sec. 
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Fig.l. Wave profile with wave period of 8 sec., in one 
wave length 
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Fig.2. Wave profile with wave period of 8 sec., in3 wave 
lengths. 

7.2. Comparison with Wiegel equation 
Using data from an observation, Wiegel (1949-1964) 
formulated relation between wave period T and maximum 
wave height H max in a wave period, i.e. 

T wieg = 15.6 .(58) 


Table. 5: Comparison with Wiegel equation 



V = 

3.0 

y = 2.97102 

T 

y z = 1.63164 

y z = 1.60095 


Hmax 

Twieg 

Hmax 

T\vieg 

(sec) 

(m) 

(sec) 

(m) 

(sec) 

6 

1,40943 

5,91305 

1,45118 

6 

7 

1,9184 

6,89858 

1,97523 

7,00002 

8 

2,50569 

7,88412 

2,57992 

8,00007 

9 

3,1713 

8,86969 

3,26522 

9,00007 

10 

3,91518 

9,85522 

4,03119 

10,0002 

11 

4,73746 

10,8408 

4,87774 

11,0002 

12 

5,63796 

11,8264 

5,80504 

12,0003 

13 

6,61695 

12,8121 

6,81286 

13,0004 

14 

7,67409 

13,7976 

7,90155 

14,0006 

15 

8,80954 

14,7831 

9,07066 

15,0006 


The comparison was done by calculating T W leg in (58) 
usingthe wave height which is the result of a calculation 


using the model, where the input in the model is wave 
period T and wave amplitude calculated using (52), so that 
the wave height that is obtained is the wave height 
maximum// maz in the related wave period. 

Table (5) shows that fory = 3, the obtained T wieg is 
almost similar with the 7'that is a wave period to 
calculate H max with the model. Whereas in y = 2.97102, 
it can be said that the obtained T wieg is equal with T. The 
result of this calculation concludes that the values of y, 
y z and equations formulated in this research are in line 
with the result of Wiegel research (1949-1964) which is 
the result of an observation. 

VIII. CONCLUSION 

If the characteristic of ideal fluid i.e. irrotational flow is 
done at Euler momentum equation, and the velocity 
potential as the product of Laplace equation solution is 
substituted, the hydrodynamic force or convective 
equation in the horizontal direction becomes zero. This 
problem can be solved using weighted total acceleration 
where there is weighting coefficient at the differential 
term against vertical-z axis and the resulted model 
producewave height that corresponds to Wiegel equation. 
Another finding that should be noticed is that the value of 
wave height is not twice the value of wave amplitude. 

A further research needed is formulating shoaling and 
breaking model by doing the weighted total acceleration 
equation, because there are many researches result in the 
laboratory on breaker height that are stated in the form of 
Breaker Height Index equation, so that the shoaling- 
breaking model and its various basic theories can easily 
be calibrated. 
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